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The unpolarized spin-flip isoscalar generalized parton distribution function (iJ" + E'^)(x,^,t) is 
studied in the large-Ai'c limit at a low normalization point in the framework of the chiral quark- 
. . . soliton model. This is the first study of generalized parton distribution functions in this model, which 

lO , appear only at the subleading order in the large- A^c limit. Particular emphasis is put therefore on the 

■ demonstration of the theoretical consistency of the approach. The forward limit of (i?" (a;, 5, 

' of which only the first moment - the anomalous isoscalar magnetic moment of the nucleon - is 

(N ■ known phenomenologically, is computed numerically. Observables sensitive to {E^ + E''){x,^,t) are 

. discussed. 
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I. INTRODUCTION 
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The study of generalized parton distribution functions (GPDs) [1] (see [2-5] for reviews) promises numerous new 
insights in the internal structure of the nucleon [6-8] . GPDs can be accessed in a variety of hard exclusive processes 
[9] on which recently first data became available [10]. Current efforts [11] to understand and interpret the data rely, 
I ' at the early stage of art, on modeling Ansatze for GPDs - which have to comply with the severe general constraints 
(«~») imposed by the polynomiality and positivity properties of GPDs [12, 13]. Model calculations [14-21] can provide 
important guidelines for such Ansatze - in particular in the case of those GPDs for which not even the forward limit 
is known from deeply inelastic scattering experiments. In this work we study the unpolarized GPD E'^(x,^,t) in the 
chiral quark soliton model (xQSM) [22, 23]. 

The model describes the nucleon in a field theoretic framework in the limit of a large number of colours Nc as a 
1^ ■ chiral soliton of a static background pion field. Numerous nucleonic properties - among others form factors [24, 25] 
as well as quark and antiquark distribution functions [26-34] - have been described in this model without adjustable 
parameters, typically to within an accuracy of (10 — 30)%. The field-theoretical character is a crucial feature of the 
xQSM which is responsible for the wide range of applicability and which guarantees the theoretical consistency of 
the approach. In Refs. [15-17] it was demonstrated that the xQSM consistently describes those GPDs, which are of 
' leading order in the large Ai'c-limit. 

In this work we extend the formalism of Ref. [15] to the description of unpolarized GPDs which appear only at 
subleading order in the large- A'^c expansion. We pay particular attention to the demonstration of the consistency 
of the approach. After checking explicitly that the model expressions for GPDs satisfy polynomiality, relevant sum 
rules, etc. we focus on the flavour combination of E°-{x, t) subleading in large- A/'c, namely the flavour singlet, which 
- being related to the spin sum rule of the nucleon - is of particular phenomenological interest. The corresponding 
flavour structure (i?" — E'^){x,£^,t), which is leading in Nc, was already studied in Ref. [15] (cf. [4]). An interesting 
issue we shall finally address is, which hard exclusive reactions are particularly sensitive to the GPD {E'" + E'^){x, ^, t). 
We observe that target single transverse spin asymmetries in two pion production are most promising in this respect. 

This note is organized as follows. After a general discussion of the properties of GPDs in Sec. II and a short 
introduction to the xQSM in Sec. Ill, we review in Sec. IV the formalism of the xQSM for GPDs of leading order 
in Nc- In Sec. V we generalize the approach to the case of GPDs subleading in Nc, and check its consistency. In 
Sec. VI we apply the approach to the numerical calculation of the forward limit of {E"^ + E'^){x,£,,t) and discuss 
phenomenological implications from our result in Sec. VII. Finally, Sec. VIII contains a summary and conclusions. 
Technical details on the calculations can be found in the Appendices. 
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II. THE UNPOLARIZED GENERALIZED PARTON DISTRIBUTION FUNCTIONS 

The unpolarized quark GPDs in the proton are defined as 

"dA 



/ 



= HO{x,^,t) U{P',s')riU{P,s) 

+E%x,U) c7(P',5')^^55^C/(P,a), (1) 

where [2:1,22] denotes the gaugc-hnk, and the renormahzation scale dependence is not indicated for brevity. The 
hght-hke vector satisfies n{P' + P) = 2. The skewedness parameter ^, the four-momentum transfer and the 
Mandelstam variable t are defined as = (P' — P)^, nA = —2^, t = A^. The antiquark distributions are given by 
H^{x,^,t) = —H'^{—x,^,t) and similarly for E^{x,^,t). The unpolarized GPDs are normahzed to the corresponding 
elastic (Dirac- and Pauli-) form factors 

1 1 
I dx H<'{x, t) = F^{t) , j dx EO{x, t) = F^{t) . (2) 
-1 -1 

The relations (2) arc special cases of the polynomiality property [2] which follows from hcrmiticity, parity, time reversal 
and Lorentz invariance, and implies that the 7V"^ Mellin moment of an unpolarized GPDs is a polynomial in even 
powers of ^ of degree less than or equal to N 



1 



da.x^-^ff'(a=,C,0)=/.2W+;,,W^2+ V e forTVeven 

° ^ 1 for odd, 



-1 

Jdx EHx, 0) = el +el^^k^ + ... + { jl, IZ Z (4) 



-JV- 



_Y^^-^ for N odd. 



As a consequence of the spin i nature of the nucleon the coefficients in front of the highest power in ^ for even 
moments N are related to each other by 

1 

/,«.W = -e«,W = |d..^-^D<'(.) , (5) 



and arise from the so-called D-tcrm _D'(z) with z = x/^ which has finite support only for \x\ < \^\ [35]. The D-term 
governs the asymptotics of unpolarized GPDs in the limit of a large renormahzation scale [4] and is related to the 
distribution of the pressure and shear forces acting on the partons in the nucleon [8] . 

In the large- A^c limit different flavour combinations of GPDs and of the D-term exhibit the behaviour [4] 

(ff" + H<i){x, ^, t) = N^f{N,x, t) , (i?" - H") (x, 1 1) = f{N,x, iV^C, t) , 
{E"+E^)ix,^,t) = NlfiN,x,N,U) , (E- " E'')ix,C,t) = Ni^fiN,x,N,U) , 
{D^+D''){z) ^Nlf{z) , {D--D^){z) ^N,f{z) . (6) 

The functions /(u, w, t) and f{z) (which are of order 0{N^) and which we do not distinguish for notational simplicity) 
are stable in the large- A^c limit for fixed values of u, v, z,t = 0{N^). E'^{x,(_,t) is systematically enhanced by one 
order in Nc with respect to H'^{x, ^, t) due to the nucleon mass Afjv = 0{Nc) appearing in the denominator on the right 
hand side of Eq. (1). As a consequence the flavour combination of H'^ix, ^, t) leading in Nc, namely {H^ + H''-){x, ^, t), 
and the flavour combination of E'>{x,^,t) subleading in Nc, namely + E'^){x,^,t), have the same order in Nc- 
This is natural from the point of view of the spin sum rule [6] 



f dx X {H'' + E''){x,£„t)=2j'^ 

•^-^ q=u,d,... 



(7) 



where 2J^ is the fraction of the nucleon spin due to (spin and orbital angular momentum of) quarks. The flavour 
singlets of E'^{x,S^,t) and H'>{x,^,t) enter the left-hand-side of (7) and contribute on equal footing in the large-A^c 
limit to J« = 0{N°). 
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In the forward limit of H'>{x,^,t) we recover the unpolarized parton distribution function fi{x) 

\im^H''{x,tt) = f!{x). (8) 
t -> 

The GPD E''{x,^,t) has also a well defined forward limit which, however, is not accessible in deeply inelastic lepton 
nucleon scattering. Phenomenologically only the first moment of E'^(x, 0, 0) is known which is given by the anomalous 
magnetic moment = -F' (0) with 

1 1 
j dx - E^){x, 0, 0) = K„ = 1.673 , j dx {E^ - E^){x, 0, 0) = = -2.033 . (9) 



III. THE CHIRAL QUARK-SOLITON MODEL (xQSM) 

The effective chiral relativistic field theory underlying the xQSM is given by the partition function [36, 37] 

Zeff = J VipV^VU exp|^j j d'^x ^{i^ - Mf/T^)vj , (10) 

where ij) and U = exp(iT°7r") are the SU{2) chiral quark and pion fields with U'^^ = exp(i75T"7r°), and M is the 
dynamical quark mass. The eff'ective theory (10) was derived from the instanton model of the QCD vacuum [37, 38]. 

An important small parameter in this derivation is the instanton packing fraction Pav/-Rav ^ 3 of the dilute instanton 
medium, where pav and i?av are respectively the average size and separation of instantons. The effective theory 
(10) contains the Weinberg-Gasser-Leutwyler Lagrangian and the Wess-Zumino term with correct coefficients and 
is valid for momenta below a scale set by the inverse of the average instanton size « 600 MeV. At this scale 
the dynamical quark mass M, which in general is momentum dependent, drops to zero. In numerical calculations 
it is often convenient to consider constant M and to apply an appropriate regularization method with a UV-cutoff 
of the order of magnitude of p~^. In calculations of (some) GPDs, however, it is important to consider momentum 
dependent M [15]. 

In the large- iVc limit, which allows to solve the functional integral over pion field configurations in the saddle- 
point approximation, the efi'ective theory in Eq. (10) describes the nucleon as a classical soliton of the pion field [22] 
providing a practical realization of the large- A^c picture of the nucleon [39] . Quarks are described by one-particle wave 
functions which are the solutions of the Dirac equation in the background of the static pion field 

^effjn) = En\n) , ^eff = -^V^fe + J° MW^' . (11) 

The spectrum of the effective Hamiltonian (11) consists of an upper and a lower Dirac continuum, which are distorted 

by the background field compared to the continua of the free Hamiltonian Hq = —ij^j'^dk + "•/^M, and of a discrete 
bound state level of energy -Eiev By occupying the discrete level and the states of lower continuum each by A^o quarks 
in an anti-symmetric colour state, one obtains a state with unity baryon number called soliton. The minimization of 
the soliton energy Esoi with respect to variations of the chiral field U yields the self-consistent pion field Uc, which for 
symmetry reasons has the "hedgehog" structure Uc{x) = exp[ierT P (r)] where P(r) is the soliton profile with r = |x| 
and Br = x/r. The mass of the nucleon is given by 

Mn = E^ipc] = minE,oi[C/] , E,oi[U] = nJe^^^ + ^ (£;„ _ En,)) . (12) 

^ E„<o ' reg 

The soliton energy is (for constant M) logarithmically UV-divergent and has to be regularized - as indicated in 
Eq. (12) and described in Ref. [25]. 

In order to include 1/Nc corrections one has to consider quantum fluctuations around the saddle point solution. 
Hereby the (translational and rotational) zero modes of the soliton solution are the only taken into account. In 
particular one considers time-dependent rotations of the hedgehog field, f7c(x) R{t)Uc{x)R^{t), where the collective 
coordinate R{t) is a rotation matrix in SU(2)-flavour space. The path integration over the collective coordinates can 
be solved by expanding in powers of the collective angular velocity = —iR^dtR which corresponds to an 1/Nc 
expansion. The latter is justified, since the soliton moment of inertia, 

AT, ^ {n\T-\j){j\T-\n) 

n, occ J "' 
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is large, namely / = 0{Nc), and the soliton rotation is therefore slow. (In Eq. (13) one has to sum over occupied ("occ") 
states n, i.e. over states with i?„ < Ei^^, and over non-occupied ("non") states j, i.e. over states with Ej > -Eiev) 
This procedure - which is referred to as quantization of zero modes - assigns to the soliton a definite momentum and 
spin-isospin quantum numbers [22] . 

The effective theory (10) allows to derive by means of path integral methods unambiguous model expressions for 
nucleon matrix elements of QCD quark bilinear operators sandwiched in nucleon states {N'\tp{zi)[zi, Z2]T'(Ij{z2)\N) , 
where T is some matrix in Dirac- and flavour space. For local observables such as, e.g., form factors the gauge link 
reduces to a unity matrix in colour space. When evaluating non-local operators in the model - as they appear in 
GPDs it is crucial that the effects of gluonic degrees of freedom, which are intrinsic in the gauge link, appear strongly 
suppressed with respect to quark degrees of freedom. This is true for twist-2 [40] (and certain twist-3 operators [41]) 
and guarantees the colour gauge invariance of the model calculation. 

If in QCD {N'\'ip{zi)[zi, Z2]rV-'(2:2)|'^) is scale dependent then the model result corresponds to low a scale of « 
600 MeV. In this way static nucleonic observables [24, 25], twist-2 quark and anti-quark distribution functions [26-34] 
have been computed in the xQSM and found to agree to within (10 — 30)% with experimental data or phenomenological 
parameterizations. In [15, 16] the approach was generalized to describe GPDs. The results of the xQSM respect all 
general counting rules of the large- iVg phenomenology. 



IV. GPDS IN LEADING ORDER OF LARGE 



It is convenient to treat the cases of flavour singlet and non-singlet quantities in the model separately for symmetry 
reasons. Let us introduce the notation 



s' s 



dX 



2tt 



(14) 
(15) 



The A^^f] are 2 x 2-matrices in spin indices which have the following behaviour in the large- TVc limit 

tr{7Vl(^=o)} = OiN^) , tr{a"'M^^=°^} = 0{N,) , 

tr{A^(-f=i)} = 0{N^) , tr{a'"A^('f=i)} = 0{NI) . (16) 

These relations do not follow from the dynamics of the model, but are group theoretical consequences of the spin- 
flavour structure of the (hedgehog) soliton field. In this sense the relations (16) are model-independent large- A^c results 
of QCD [42] which are consequently respected in the xQSM. 

In order to evaluate the bi (nucleon- )spinor expression on the right-hand-side of Eqs. (1) one has to consider that 

in the large Nc limit the nucleon is heavy, Afjv = 0{Nc), and the kinematics becomes non-relativistic. In particular 
we have for the components of the momentum transfer A* = 0{N^) and A° = 0{N~^). Thus the hierarchy 



holds Mn » |A'| > |A°|, while t 



0{N^) and ^ = -A^/{2Mn) = 0{N-^). Here we have chosen 



n'^ = (1, 0, 0,— 1)/Mjv. Evaluating consequently the right-hand-side of Eq. (1) in this large-A?'c kinematics yields for 

the flavour-singlet and non-singlet case, respectively. 



(1=0) 



^3kl 



2 5.,.,' 



t 



N 



{H^ + H'') {x, t) + {E^ + E''){x,^, t) 



{H- - H''){x,^,t) + -^{E^ - E'){x,£,,t) 



N 



■ia'),,,{E"-E^)ix,^,t) 



From Eqs. (17, 18) simple relations follow for the leading (in the large- iVc counting) GPDs 



{H^ + H''){x,^,t) = ^tr{M^'-°^}, 



iE^-E'')ix,^,t) 



iM, 



N€' 



-Sbm 



2 A 



tr{CT"M(^=i)} 



(17) 

(18) 

(19) 

(20) 



where "tr" denotes the trace over spin indices and A^^ = 
however, appear only in combination with the leading ones 



2 /A3^2 



(A3)2 



t — 4:M^(^'^. The subleading structures, 



{H'' + H'^){x,^,t) + {E'' + E'^){x,^,t) = 2 tr{a'"A^(^=°)} 

(ff" - H''){x,^,t) + ^{E-- E'^){x,i,t) = ltr{M<-'-'^ . 



(21) 

(22) 
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The structures in the square brackets of Eqs. (17, 18) and on the left-hand-sides of Eqs. (21, 22) are analogs of 
the relations between the electric and magnetic form factors G\.(t), G\^{t) and the Dirac- and Pauli-form factors 
F^{t), F^{t), which are given by 

GUt)=FUt)+FUt), Gl{t)=FUt) + ^F^{t). (23) 

The model expressions for A^(^='^'i) in leading order of the large- iVc limit were derived in Ref. [15] and the following 
results for (ff" + H'^){x, i, t) and (£;" - E"^) were obtained using Eqs. (17, 18) 

" 7l,OCC 



Jz°{xMn-E„) 



X (n| (1 + 7°7^) exp{-i^f ) exp(iAX) exp(-z^p^) \n) , 



(24) 



3(A-^)2 



2tt ^ 



X (n| (1 + 7°7^)(t X A)^ exp(-i^p^) exp(iAX) exp(-i^p^) |n) . 



(25) 



There are equivalent expressions with ^Yliocc ^ ^Snon' where the summation goes over non-occupied states 
En > Eicv The possibility of computing in the model quantities in these two independent ways is deeply related 
to the locality properties of the model [26]. Numerical results for (24, 25) were presented in Ref. [15], and further 
discussed and reviewed in Ref. [4] . 

Let us emphasize that according to Eq. (8) in the forward limit the right-hand-side of (24) reduces to the model 
expression for (/f + ff){x) [15]. (7?" -I- H'^){x,S^,t) and {E"^ — E'^){x,^,t) are correctly normalized to the respective 
form factors [15], cf. Eq. (2), and they satisfy the polynomiality conditions in Eqs. (3, 4) [17]. 

It is worthwhile mentioning that the coefficient of the highest power in ^ of even Mellin moments of (S" — £''') {x, ^, t) 
in leading order of large Nc is zero in the xQSM. Also this result is not a dynamical feature of the model, but rather 
a group theoretical consequence of the soliton symmetries which ensure the correct large- A^c counting for the flavour- 
nonsinglet D-term. In fact, if this coefficient were not zero, then (£"-£'^) = 0{N^) would imply {D'^-D'^) = 0{N^) 
in conflict with the counting rule which states (D" - Z)**) = 0{Nc), cf. Eq. (6). 



GPDS IN SUBLEADING ORDER OF LARGE Nc 



One has to consider 1/A^c (rotational) corrections to Al^^^"-^) in order to study the subleading flavour combinations 
(iJ" — H'^)(x,^,t) and {E^ + E''-){x,£_,t). In order to simplify the notation let us introduce 

(26) 
(27) 



EM{x,^,t) = {H^+H^){x,i,t) + {E" + E^){x,^,t) ~ 0{N!) , 
He{x, C, t) = {H^ - H''){x, ^, t) + {E- - E'')ix, U)^- 0{N,) . 



Here the respective already known (see above) flavour combinations of H'^{x, ^, t) and E'^ix, ^, t) leading in Nc appear. 
Thus it is sufficient to focus on the "new" objects EM{x,^,t) and HE{x,^,t). 
The model expressions for HE{x,^,t) and EM{x,(,,t) (for a proton) read 



HE{x,^,t) = - 



MnNc 
12/ 



'f L VTn,oc< 



-E 



-iz^'E, 



,occ 



1 



+ ■ 



d 



Em - Ej OxMn 



EM{x,£,,t) 



E 

Tn,occ 

7 OCC Ji^^ll 

X (mlr^lj) (il r"(l + 7V) exp(-i^V/2) exp(iAX) exp{--iz°f 12) \m) , 

1 d 



-iz"E„ 



Jz°xMm 



(28) 



21 



dzO 



J2 e-"°^'"-J2 



m.occ 
J, all 



m,all 

J, OCC 



X (mlrli) 01 (1 + 7%") exp(-i^V/2) 



Em - Ei OxMn 



E 

m.occ 
3, all. 



-lz"Err 



^iz xMn 



exp(iAX) exp(-i^V/2) |m) . 



(29) 



There are equivalent expressions with opposite sign where the summations go over non-occupied states. Considering 
the generalizations due to the oS'-forward kinematics [15, 16], the derivations of Eqs. (28, 29) closely follow the 
derivations of the model expression for the flavour non-singlet unpolarized distribution (/" — fi){x) [30] and the 
flavour-singlet helicity (5" + gi){x) and transversity (/i" + hf){x) distribution functions [33]. In the following we check 
explicitly the theoretical consistency of the expressions (28, 29). 
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A. Form factors and polynomiality 



Let us first verify the correct normalization of HE{x,^,t) and EM{x,^,t). In order to integrate Eqs. (28, 29) over 
X it is convenient to substitute x —t y = xMj^ and to extend the y-intcgration range [—Mj^,Mn] to [—00,00] in 
the large- A^c limit. The derivative in x drops out and after cancellations in the curly brackets in Eqs. (28, 29) only 
summations over occupied states |m) and non-occupied states \j) and vice versa remain, which can be combined by 
exploring model symmetries, see App. A. We obtain 

1 

/ dxHE{x,^,t)=-^Yl ;^^Hr1j)0|r"exp(zAX)|m)^Gg=^)(i), 

m.occ J 



/AT A/f j^abcAa 
dxEM{x,U) = E ;g;;^Hr''|i)(i|7Vexp(iAX)|m) ^3Gi^=°)(t) , (30) 

where we identify the model expressions for the electric isovector G^~^^(t) and magnetic isoscalar G^lf^\t) form 
factors [24, 25]. [The factor 3 in the second line of (30) appears because for GPDs the notion of (non-) singlet commonly 

refers to quark flavours, e.g. (i?" + H'^). In contrast, in the case of form factors it refers to proton and neutron, e.g. 

G^"°^ = + GIj.\ From (23, 26, 27, 30) and recalling that the leading large- A^c GPDs (il" + H'^){x,^,t) and 
(£^" — E'^){x,^,t) are correctly normalized [15], we find for the subleading GPDs in agreement with Eq. (2) 

1 1 
I dx (ff" - H'')ix, 1 1) = {F^ - Ff){t) , J dx {E^ + E''){x, ^, t) = {F^ + Fi){t) . (31) 

-1 -1 

In App. B it is explicitly demonstrated that the higher moments in x of He{x, ^, t) and Em{x, ^, t) are even polynomials 
in ^ according to (3, 4). In fact, we find that in even moments the coefficients in front of the highest power in ^ have 
opposite sign in (£" -|- E'^){x, ^, t) and (ff" -|- H'^){x, ^, t) in accordance with the relation (5). 

B. Spin sum rule 

For the second moment of Em{x, ^,t) at t = we obtain in the xQSM, cf. App. C 

da; xEm {x,^,0) = 2S'^ + 2L'^ = 2S^ = 1 , (32) 



/: 



where = ^ is the total spin of the nucleon, and and are the respective contributions of the spin and orbital 
angular momentum of quarks and antiquarks to the nucleon spin. The result in Eq. (32) is consistent. In the effective 
theory the contribution of quark and antiquark degrees of freedom, J'3 = 5"^ -|- L*^, must account entirely for the 
nucleon spin, since there are no gluons in the model. Thus, the result in (32) means that the spin sum rule (7) is 
consistently fulfilled in the model: 



1 

jdx xiH"" + + E"" + E^){x,^,0) = 2J'^ = 1 



(33) 



We note that in the model 25*^ = g'^^ = 0.35 where 17^ denotes the isosinglet axial coupling constant [33]. Thus, in 
the xQSM 35% of the nucleon spin are due to the quark spin, while 65% are due to the orbital angular momentum 
of quarks and antiquarks. 

In a gauge theory it is not possible to separate unambiguously spin and orbital momentum [43, 44] . It also is by no 
means clear that the identification of the model expression for quark orbital momentum in Eq. (32) is unambiguous. 
It is interesting to note that one arrives at the same result by intuitively identifying = e^^^x^p^ with the orbital 
angular momentum operator of the effective theory (10), cf. [45]. However, one has to keep in mind that this issue c;an 
rigorously be studied in the instanton vacuum model from which the xQSM was derived. Noteworthy in this context 
is that in the chiral quark soliton model arises from the realization of the axial anomaly in the instanton vacuum 
[40]. 

Since in the chiral quark-soliton model also the total momentum of the nucleon is carried by quarks and antiquarks 
only [26], i.e. = 1, we obtain for the second moment of (£■" -|- E'^){x, 0, 0) the result 

1 

jdx a;(£;" + E'^){x,Q,{)) = {2J^ - M^) = . (34) 
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Interestingly, the result (34) holds also in QCD in the asymptotic limit of a large normalization scale /U — > oo. This 
happens because M'^ and 2J^ have the same asymptotics [43] 

lim MQ = lim (2J«) = (35) 

where Np is the number of quark flavours. It should be stressed, however, that at a low scale /i ^ 600 MeV in the 
model Af^ = 27'^ = 1 is far from the asymptotic values in Eq. (35). 

Eq. (34) is the general prediction of any model lacking explicit gluon degrees of freedom, which is able to consistently 
describe the nucleon in terms of quark and antiquark degrees of freedom - such as the xQSM. One can show, using 
nic;thods of theory of the instanton vacuum [40] that the gluon contribution to the nucleon momentum and angular 
momentum is parametrically suppressed by the packing fraction of the instantons in the vacuum. Therefore in order 
to obtain non-zero gluon contributions one has to extend xQSM beyond the leading order in the instanton packing 
fraction. 

Noteworthy, it has been argued that momentum and angular momentum should be equally distributed among 
quarks and gluons at any scale, not only in the asymptotic limit (35), an observation which can be reformulated as 
the absence of an anomalous gravitomagnetic moment of the nucleon [46] . 



C. Forward limit 



Finally let us discuss the forward limit. From (27) we see that He{x, 0,0) = (H'^ — H'^){x,0,0). Indeed, taking 
A in Eq. (28) and making use of the hedgehog symmetry we recover, in agreement with Eq. (8), the model 
expression for the flavour non-singlet unpolarized distribution function [30] 



He{x,o,o) = -nri- 2^ 



12/ 



d 



3, all 



Ej — E„ 



dxM, 



N 



{mn){j\T''{l + lW)KxMN-Em-f)\m) = (/r-/f)(x). (36) 



The forward limit of Em{x, t) {E'"' + E'^'){x, 0, 0) + (/" -|- ff)(x) contains a contribution which is a priori not 
known, namely {E^ +E'^){x, 0, 0). Therefore we derive here the model expression for Em{x, 0, 0) which we shall study 
below in detail. The limit ^ ^ is regular^, however, to complete the forward limit we have to consider with care the 
limit A I — > of the structure 



fZjk \3 A J \3 Am 

-^2^ exp(iAiXi) = e^^'^ + i^^^'X^ + 0(A'=) . 



(37) 



The flrst term in the above expansion (for small but non-zero A_l) yields a vanishing result when inserted into 
Eq. (29) due to the hedgehog symmetry. The second term from the expansion in (37) yields the only contribution 
which survives the forward limit in Eq. (29). Making use of 



/\3 Am 1 1 



(38) 



we obtain upon use of hedgehog symmetry 



Em{x, 0,0) 



41 



E 

m.occ 
J .all 



Ej — E„ 



dxM, 



N 



(mlr'^li) 01 (1 + 7V) S{xMn - S„ - p') |m) . (39) 



As a final check we integrate Em{x, 0, 0) in Eq. (39) over x and recover the model expression for the isoscalar magnetic 
moment [49] 



dxEM{x,0,0)= 2^ 



1 



Ej — Er, 



|r^li)0l7°(7xX)3|m) = 3M(^=°^ 



(40) 



Of course, by limiting oneself to § = one drops important physics. E.g., the skewedness parameter is related to the correlation length 
between partons, see [47] and references therein. The present calculation unfortunately cannot shed any light in this respect. 
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Imco 




FIG. 1: The contour C of the w-integral in Eq. (41). 



VI. THE FORWARD LIMIT OF (£" + E'^){x, 5, t) 

In order to compute the forward limit (£" + £:'')(.t, 0, 0) = Em (x, 0,0) - (/f + we have to numerically 

evaluate the model expression (39) for Em (x, 0,0). The distribution function (fi + ft){x) was already studied in 
[26, 28, 29]. The numerical methods needed for that were developed in Refs. [28-30] in the context of usual parton 
distributions functions in leading and subleading order of the large- iVc expansion. We restrict ourselves to the forward 
limit since this numerical technic does not allow a full calculation of GPDs for all values of $, and t. 



A. Gradient expansion and chiral enhancement of GPDs 



Before studying the model expression (29) numerically in the model we shall consider the gradient expansion of 
EM{x,£,,t). This expansion consists in expanding model expressions in powers of the gradients of the (static) chiral 
field VC/. Such an expansion would quickly converge if the soliton field were slowly varying, i.e. if the soliton were 
large VU ~ l/i?soi M where i?soi is the scale characterizing the soliton size. However, for the physical soliton 
solution Rso\ 1/M. Nevertheless the gradient expansion is instructive and allows - among others - to study the 
UV-behaviour of the model expressions. Also as we shall see below the gradient expansion allows to obtain strong 
enhancement of GPD Em{x, ^, t) in the region of small x. This enhancement is related to the effect of the pion cloud. 

The model expression for EM{x,^,t) can be rewritten equivalently in the following way 



EM{x,^,t) = 



f dcjd^p ^ 



27 Ai 



(5(u;+p^-a;MAr)(l + 7V) 



u-H Lo-H 



+ 



dxM. 



N 



5(uj + p'^ - xMn){1 + 7 7 )t 



H 



(41) 



where "tr-^^i?" denotes the trace over Dirac- and flavour indices. The contour C is defined in Fig. 1. Closing the contour 
C in the upper half of the complex cj-planc yields the expression in Eq. (29), closing it in the lower half plane one 
obtains the equivalent expression whore the summation goes over non-occupied states m, cf. the sequence of Eq. (29). 

Expanding the expression in Eq. (41) in powers of V(7 one obtains^ already from the zeroth order in (VJ7) the 
following result for the continuum contribution 



ET\x,^,t) = 



1 



3/(k+f)/(k-f) 



47Aiy (2.).-^"' ,.-A^.3^^3(^,^ 
[ ejxMN + k^)Q{^ - x) - Q{-xMn - k^)e{x - Q 

\ fc3 - ^Mn 



2 For that we rewrite the operator l/(aj - H) = i(u) -|- H)l{u? - H'^)+ ^(u)'^ - H^)-^{ui + H) where H'^ = + M'^ + iM'Y^U~"i , and 
expand it in powers of iMf^V^^ . Such a "symmetric expansion" of l/(a; — H) ensures the hermiticity of the operator also in the case 
when the series in Vi7 is truncated. 



9 



2^Mjv J 
where /(q) = /(|q|) is one of the functions describing the Fourier transform of the U{x)-&eld 

d^x ([/(x) - 1) e-'^- = 5(q) + ^/(q) ; /(q) ^ 4^ / dr r2ji(r|q|) sinP(r) . (43) 

111 ^0 

Note that the property EM{x,^,t) — Em{x, —S,,t) holds manifestly for the result in Eq. (42). 

In the derivation of Eq. (42) we encountered a logarithmic UV- divergence proportional to M^. This UV-divergence 
can be regularized, e.g., by means of a Pauli-Villars subtraction as 

i;^"*(x,$,i)reg = E'jS'^\x,C,t;M) - i;-"t(a;,e,t;Mpv). (44) 

We removed the dependence of the result on cutoff and regularization scheme in favour of a physical parameter, 
namely the pion decay constant f„ = 93 MeV given in the effective theory (10) by the (Euclidean loop) integral 



J {2^[ 



)4 (p| + M2)2 



(45) 



which exhibits a similar UV-behaviour and can be regularized analogously. In the first order of the gradient expansion, 
indicated only symbolically as 0(VU) in Eq. (42), we also encounter a logarithmic divergence which can be removed by 
means of (44). Still higher orders in the gradient expansion yield finite results - as one can conclude from dimensional 
counting. (Or from the fact that the isoscalar magnetic moment appears at the order (V?7)^ and is UV-finite, cf. 
Ref. [23].) 

Eq. (42) may give at best only a rough estimate for EM{x,S,,t). Nevertheless it contains already main features of 
the total result. Apart from the UV-properties, which we discussed above, the small x-behaviour in the chiral limit 
is of interest since the exact numerical calculation is performed under such conditions. The soliton profile behaves at 
large distances as 

B 3 a'-^^ 

Pir) = — o (1 + m^r) exp(-m^r) , B = for r » \/l (46) 

(3) 

with B related to axial isovector coupling constant g\' = 1.26 by means of equations of motion. The large distance 
behaviour (46) of the soliton profile translates into the infrared-behaviour of the function /(q) as 

/(q) = -^^j^ for small q2,m^«B-^ (47) 

Note that B~^ < {An Jt^)"^ ~ IGeV^ which is a typical scale for chiral symmetry breaking effects. Apparently the 
limits m„ — > and q — > do not commute.^ In the forward limit in Eq. (42) this means 



£^m(x,0,0) = gj-^27r)2 J'^'i^fw " 



35 



(3) 



\xMn\ o for ^ 0. 




for = 0, 

(48) 



Note that the soliton moment of inertia could be eliminated in favour of the A-nucleon mass-splitting as Ma — Mm = 
3/(2/). 

Although derived in the leading order of the gradient expansion of model expressions Eq. (48) is of more general 
nature as it is related to the effect of pion cloud. Therefore it is should be possible to derive the above expression by 
means of chiral perturbation theory. Another interesting question is to find phenomenological manifestations of the 
chiral enhancement. This would open new exciting possibilities to investigate chiral symmetry breaking in the hard 
exclusive processes. 



The iioii-coiiiiimtativity of such hiuits is known, e.g., from studies of the polarized photon structure function, where the virtuality of the 
photon plays the role of the momentum q. In that case it is known that the correct order of limits is to take first the photon virtuality 
to zero, and only then to consider m„ [50]. However, in our context it is not clear in which order the limits should be taken, though 
phenomenology may suggest first to keep mn finite. 
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Numerical calculation 



As any quantity in the model, Em{x, 0, 0) is composed of a contribution of the discrete level and of the continuum 
contribution defined, respectively, as 



MlNr. 



AI 

MlNr. 



Ej ,all 
j ^ lev 



d 



Em (a;,0,Oj = — 77— e 2^ 



Ej-E\ev BxMn 

2 d 



41 

MlN, 



47 



Em>0 
Ej ,all 



Ej - Em dxMN 



d 



Ei - Em dxM, 



N 



(levlr'^lj) (j| (1 + 7°7') S{xMi, - E^ - f) |lev) 
(m|r'=|j) (jI (1 + 7 V) 5{xMn -Em- p') K |m) 

(mlr'^li) 01 (1 + 7 V) S{xMm -Em- f) \m) 



(49) 



For a constant (momentum-independent) Af the continuum contribution to Em{x, £,,t) is UV-divergent. It can be 
regularized by means of a single Pauli-Villars according to Eq. (44). In the context of parton distribution functions 
the Pauli-Villars method is the preferred regularization scheme because it preserves fundamental properties of parton 
distributions such as sum rules, positivity, etc. [26]. It should be noted that the contribution of the discrete level, 
which is always finite, must not be regularized [29]. 

The numerical method to evaluate £^M(a;,0,0) consists in placing the soliton in a large but finite spherical box, 
introducing free basis states (cigcnstatcs of the free Hamiltonian, sec Eq. (11) and below), and discrctizing the basis 
by imposing appropriate boundary conditions following Kahana and Ripka [48]. The full Hamiltonian (11) can be 
diagonalized numerically in this basis [49]. With the eigenfunctions and eigenstates, $n(x) and En, one is then in a 
position to evaluate Eq. (49). This is most conveniently done by converting the model expression (49) into a spherically 
symmetric form [28] . For that we introduce in (49) a unit vector a such that the original expression is recovered for 
a = e(3), 



e3jfe(i + ^0^-) ^(^^^ -Em- f) Xi^{l+ 7°a • 7) S{xMn - ^™ - a • p) (a x X)* 



(50) 



Then we average over all possible orientations of the vector a. This yields 



^M(a;,0,0) = - 



47 



m.occ 
m^3 



Em - Ej OxMn 



{m\T''\j){j\ 



(pxX) 



k h{\p\,Em,x) 



+(7°7 X X) 



where Ii{\p\) = 



fc (7o-72)(lpl,j^^,a:) , ,,^^,fc, 0,. (372 -7o)(lpl,£;„,a:) 



2|P| 

{xMn - Em Y 



(pxX)'=(7^.p). 



2|P|' 



\m} , 



21p 



Ml 



e(lp] - \xMn - E„ 



(51) 



Although non-commuting operators appear in (51) nevertheless the final result is a hermitian operator. The averaging 
procedure must, of course, preserve hermiticity. It is this step which prevents us from tackling also the problem to 
compute numerically Em{x, £,,t) for non-zero ^ and t. For finite A the appearance of an additional direction, namely 
along Aj^ perpendicular to the 3-direction, makes the resulting expressions unsuitable for a numerical evaluation. 

Since we work with a discrete basis, the ©-functions in (51) would lead to discontinuous functions. Therefore we 
convolute the true (but due to finite size effects discontinuous) function with a narrow Gaussian 



00 

7;^(x,0,0) = -^ y dx' e-(^-^')'/''' 7i;m(x',0,0) 



(52) 



The width 7 has to be chosen adequately according to the spacing of the discrete basis states [28]. The effect of 
smearing can be removed by a deconvolution at the end of the calculation. 

Fig. 2a shows the level and continuum contribution to Em {x, 0,0) computed with M ~ 350 MeV in a box of the 
size 7)box = 18 fm with 7 = 0.08. The continuum contribution is evaluated in the two equivalent ways, Eq. (49). This 
explicit check of the equivalence of summations over respectively occupied and non-occupied states is an important 
check of the numerical method, and provides a measure for the numerical accuracy. As can be seen from Fig. 2a the 
two calculations yield the same result within few percent. 
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FIG. 2: Em (a;, 0,0) = (_B" + E'' + + //'*)(:r, 0, 0) as function of x from the xQSM at a low scale of about ^l ^ 600 MeV. At negative x 
the curves show minus the antiquark distributions. Fig. 2a shows the (smeared) discrete level and continuum contribution. The latter is 
evaluated in the two equivalent ways in Eq. (49), namely by summing respectively over occupied and non-occupied states. In Fig. 2b the 
effect of the Gaussian smearing is removed. These are the final results for respectively the discrete level and the continuum contribution 
to Em{x, 0,0). 

We observe that the equivalence is spoiled by unphysical boundary effects in the chiral limit, when the soliton profile 
P(r) oc 1/r^ at large distances. In order to restore the equivalence it is necessary that the profile decays faster at large 
r, e.g., as P(r) oc cxp(— m^r)(l + m^r)/r^ (the numerical parameter m'^ is not to be confused with the physical pion 
mass discussed above) . Since the issue of computing self-consistent profiles with finite pion masses in the Pauli-Villars 
regularization is not yet solved, cf. [51] for a discussion, we use the chiral self-consistent profile computed in [29] up 
to some ta and continue it for r > with an artificial Yukawa-tail-like suppression of the above kind. We find the 
results practically independent of and m'^ in the range r^i = (4 — 8) fm and m'^ = (100 — 200) MeV. This proves 
that what matters in this context is only a sufficiently small value of the profile at the boundary r = -Dbox, and makes 
the the extrapolations ta — * oo and/or — » superfluous, which in principle would be necessary to remove any 
dependence of these numerical parameters. The equivalence demonstrated in Fig. 2a is achieved in this way. This 
procedure was applied successfully to cure an analogue problem in the calculation of the transvcrsity distribution [34] . 

Let us comment on the effect of smearing, which is negligible whenever one deals with a continuous function such 
as the contribution of the discrete level, cf. Fig. 2a and 2b. However, in the case of the continuum contribution the 
smearing "hides" an l/x singularity which becomes apparent only in the final result after the smearing is removed. 
This is done by Fourier transforming i?U-"(a;, 0, 0) of Eq. (52), dividing out the smearing Gaussian, and re-Fourier 
transforming - which yields the final results shown in Fig. 2b. 

We remark that in the xQSM the parton distribution functions (and forward limits of GPDs) do not vanish for 
> 1. Instead they decay as exp(— const iVc a; ) at |a;| > 1 [26] which is, in fact, numerically very small even for 

C. Discussion of the results 

We observe in £'M(a^, 0,0) that the discrete level contributes predominantly to the distribution of quarks rather 
than antiquarks. The continuum, however, contributes nearly equal portions to quark and antiquark distributions. As 
a result the isoscalar magnetic moment receives a negligible contribution from the continuum. This strong dominance 
of the level contribution in the isoscalar magnetic moment was observed in earlier calculations [24] . The result (in 
units of the nuclear magneton) 

1 

/x(^=°) = {iiP + /i") = i y^da; Em{x, 0, 0) = 0.65 vs. 0.88exp , (53) 
-1 

agrees with the experimental value to within 25%, i.e. to within an accuracy typical for xQSM results [25]. 

We can interpret xEm{x, 0, 0) = x(7?" + H'^ + + E'^){x, 0, 0) (and analog for antiquarks) as the distribution of 
spin in the nucleon in the following sense. This quantity tells us how much of the nuclcon spin is due to quarks (and 
antiquarks) from which region in x. In Fig. 3a we see that the contribution due to quarks is concentrated around 
X « 0.3, while the contribution of antiquarks has its maximum in the region x « 0.1. It is worthwhile stressing 
that a sizeable fraction of the nucleon spin is due to antiquarks already a the low scale of the model. We find the 
nucleon spin distributed among antiquarks and quarks as antiquarks : quarks =1:4. Integrating xEm{x,Q,0) over 



0. 
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distribution of spin distribution of momentum spin vs. momentum 




FIG. 3: Comparison of the "spin" and "momentum distribution" in the rmcleon in the xQSM at a low scale of about /i = 600 MeV. 
Fig. 3a shows the respective contributions of quarks and antiquarks to the rmcleon spin as distributions in x, defined as a;i?M(3;, 0, 0) = 
x{E^ + E'^ + H^ + H'^)(x,Q,Qi) (and analog for antiquarks). The sum of the integrals over x of the quark and antiquark contributions yields 
2J'3 = 1. Fig. .3b shows the distribution of the momentum in the nucleon, + ff){x), for quarks (and analog for antiquarks). The 

sum of the integrals over x of the quark and antiquark contributions yields iW^ = 1. In Fig.2c the total spin and momentum distributions 
are compared. 

X and summing up the contributions of quarks and antiquarks we obtain unity (within 2% numerical accuracy) - in 
agreement with the spin sum rule in the model, Eq. (33). 

It is instructive to compare - within the model - the spin distribution to the momentum distribution of quarks and 
antiquarks in the nucleon. First we note that antiquarks arc sizeable also in the momentum distribution, however, 
less pronounced than in the case of the spin distribution, see Fig. 3b. The nucleon momentum is distributed as 
antiquarks : quarks « 1 : 10. Parameterizations performed at comparably low scales [52] indicate that about 30% of 
the nucleon momentum is carried by gluons already at low scales. Thus, the distribution of momentum in the model 
seems to overestimate the contribution of the quarks and antiquarks by about 30%, which is within the expected 
accuracy of the model [31]. It will be interesting to see whether the situation is similar in the case of the spin 
distribution. 

Finally, comparing the momentum and spin distributions. Fig. 3c, wc observe that the spin distribution is somehow 
softer, i.e. shifted towards smaller x. 

In Fig. 4a, we present the final results for x{E'^ + E'^){x,Q,Q) and x{E^ + E'^){x,Q,Q). We observe that xi^E"" + 
£"^)(x, 0, 0) is positive at small x, changes sign at x ~ 0.18, and is sizeable and negative at large x, taking its maximum 
around x « 0.5. In contrast, x(i?" + E'^){x, 0, 0) appears everywhere positive and rather large at smaller values of x. 
For the second moment of (i?" + E'^){x, 0, 0) we obtain zero (within numerical accuracy) and confirm numerically the 
sum rule in Eq. (34), which - let us stress it - must be fulfilled in any model which lacks explicit gluon degrees of 
freedom but describes the nucleon consistently in terms of quark- and antiquark degrees of freedom. 

In our calculation x(i?"-|-i?'^)(x, 0, 0) results from the cancellation of two much larger quantities, namely xEm{x, 0, 0) 
and x(/" + ff){x), see Figs. 3 and 4, such that one could worry whether a small change in parameters could lead to a 
much different final result. The only parameter, which could be in principle varied^, is the constituent quark mass M. 
A calculation with M = 420 MeV yields a somehow larger result, see Fig. 4b. However, what is more important in our 
context, it fully confirms the basic features of the calculation with M = 350 MeV in Fig. 4a, which is our final result. 
Fig. 4c shows the model results LO-evolved [58] to = 5GeV^ which is a typical scale in experiments. Hereby we 
assume the gluon GPD E^ {x,0,0), which mixes with (£■" + E"^ + E^ + E'^){x,0,0) under evolution, to be zero at 
the initial scale taken to be — 0.36 GeV^. This is justified by the instanton vacuum model, where twist-2 gluon 
distributions appear suppressed with respect to the respective quark distributions. (The meaning of the dotted line 
in Fig. 4c is explained in the next section.) 



* Wc work here with the value M = 350 MeV which follows from instanton phenomenology [37, 38] and was used in calculations of parton 
distributions and GPDs [15-17, 26—33]. However, in numerous model calculations M was allowed to vary in the range (350 — 450) MeV 
[25]. The value M = 420 MeV was somehow preferred because it reproduced exax;tly the delta-nucleon mass-splitting within the 
proper-time regularization [25]. 
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FIG. 4: Fig. 4a shows the final result for x{E^ + E''-)(x, 0, 0) and a;(E" + E''-){x, 0, 0). For sake of comparison, see text, Fig. 4b shows 
the result one obtains with a constituent mass M = 420 MeV. Fig. 4c shows the final model result (Fig. 4a) LO-cvolvcd from the initial 
scale taken as /i'^ = 0.36 GeV^ to the scale = SGeV'^ typical in experiments. For comparison the Ansatz of Ref. [4] is shown for 
quarks (dotted line). The corresponding antiquark distribution was assumed to be zero in Ref. [4]. Note that the effect of sea quarks was 
simulated in [4] by a <5-function not visible in the plot. 



VII. PHENOMENOLOGICAL IMPLICATIONS 



Aziimithal asymmetries in hard exelusive produetion of vector mesons or pion pairs from a transversely polarized 
target are among the observables which are most sensitive to E'^{x,^,t) [4]. Let us define 



1 /;d(/.aW-/J"d0 c7(</.) 
\S±\ J^^dcPa{ct>) 



AMixB,t,Q^)= /° ^ X , , ,/ ^ (54) 



where a is tlie cross section for the process jKq) + P{p) M + P{p') with jl{q) denoting the deeply virtual 
longitudinally polarized photon with momentum g. P is the incoming (outgoing) proton with momentum p{p'), and 
M denotes the produced longitudinally polarized vector meson or the pion pair. For notational simplicity we omit 
to indicate that a is differential in t = (P — P')"^, = — g^, xb = Q'^/{2Pq) which is related to the skewedness 
parameter as xb = 2^/(1 + ^) in the limit Q"^ ^ oo, and the angle </> between the transverse proton spin S± and the 
plane spanned by the virtual photon and the produced meson. 

For not too small Xb the gluon exchange mechanism can be neglected in a first approximation and the asymmetry 
Am reads [4, 53] 



2|Ai| lm(B*C) 

i-M ~ 



nMM W{1 e) \C\-'{e + mT)- Re{BC*)2e 
where B and C are given for the respective final state M by 



(55) 



Bpo = [\x{euH^-edH'^)( 1— + — ^ 

J-i \x-(,+ie x + £, 

Cpo = [ dx {euE^ - edE'^) ( 1— + —1- 

J-l \X-^ + lSX + ^' 



le 



1 ^11 



= / dx {euH^ + CdH'') { . 

' ^^x — t^+iex + t,—i£ 



j dxieuE"" 



a = / dx {euE-^ + BdE") 



1 1 



B^o^o = j dx{euH'' + edH'^)(^-^ 



X — S,-\-i£ a; + ^ — 
1 



X — x + 



C^o^o = [\x{euE^ + edE'')( 1— ] . ] 



(56) 



An advantage of choosing such target spin asymmetries to access GPDs is that the (poorly known) meson distribution 
amplitudes cancel out. 
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As we computed here only the forward hmit and are particularly interested to see the impact of our result in 
comparison to what previously was discussed in literature, we shall use the method of Ref. [54] to generate ^- and 
f-dependence from a given forward limit according to 

E\x, ^, t) = EIj,{x, - D-term (57) 

with the D-term as predicted in the model [4], F2{t) from [55, 56] and the so-called double distribution given by 

/•I /-i-l/?! 

EIj^{x,£)= dp daS{x- (3- a^)Ei{l3, 0,0) h{j3, a) . (58) 

J-i J-i+\0\ 

The Ansatz (57, 58) satisfies the polynomiality condition (4) for any "profile fimction" h{(3,a), e.g., for 

r(2b+2) [a-mr-a']'' 

~ 22b+ir2(6 +1) (1 - 1/31)26+1 ^^^> 

where wc shall choose 6 = 1 for our analysis, cf. Ref. [4]. For (_£" + E'^){x,0,0) we use the result obtained in this 
work, LO-evolved^ [58] to Q2 = 5 GeV^ cf. Fig. 4c. For - E'^){x,0,0) we shall employ the Ansatz 

(£« - E'')ix, 0, 0) = Ci(A".,, - /tai)(a;) + c^Six) (60) 

inspired by the model results for (E" — E'^)(x,0,0) [4]. In Eq. (60) ff^^i denotes the respective valence quark 
distribution. The parameters ci,2 are fixed from — k"^ = 3.706 and, in order to discuss everything in terms of results 
from the model, J"" - ^ 0.2 found in the xQSM [4]. 

We shall compare the results obtained with this Ansatz, which is consistently based on predictions from the xQSM 
to the results presented in Ref. [4] where, in lack of better knowledge, the Ansatz (60) was assumed to hold also in the 
flavour singlet case. Such an assumption is, however, not supported by our direct calculation as it is demonstrated in 
Fig. 4c. 

The GPD H'^{x,^,t) is modeled in an analog way on the basis of the parameterization [59] for (x) as described 
in Ref. [4]. 

As can be seen from Fig. 5 the asymmetry has little sensitivity to the Ansatz used for {E"^ + E''-){x,0,0) in the 
case of u>^ or p° production. We find that the process most sensitive to (S" -|- S'^)(x,0,0) is the production of a 7r° 
pair. The effect of different Ansatze can cleanly be distinguished. Experiments at HERMES or JLAB could provide 
interesting insights. However, for an unambiguous quantitative analysis of the data it will be necessary to consider 
systematically power and NLO corrections. 



® One could also first construct the GPD at the low scale in the above describe way, and then evolve it to the relevant scale. The difference 
due to the non-coinmutativity of these steps is, however, small [57] — fax smaller than other uncertainties in our model calculation as 
well as corrections due to the neglect of NLO-corrections, possible power corrections, etc. We shall neglect it here. 
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VIII. SUMMARY AND CONCLUSIONS 

We have presented in the framework of the flavour-SU(2) version of the chiral quark-soUton model a study of 
unpolarized GPDs, which appear at sublcading order in the large- iV^ limit, namely the flavour combinations {E" + 
E''-){x,^,t) and (if" — H'^){x,^,t). For that we generalized the methods of Ref. [15] developed to study the leading 
large- A^c GPDs, i.e. the respectively opposite flavour combinations. 

Interestingly, the GPDs (iJ" + H'^)(x, ^, t) and {E" + E'^){x, ^, t) are of the same order in the large- iVc counting but 
handled differently in the model. The former appears already in leading order, while for the latter one must invoke 
1/Nc (rotational) corrections to obtain a non- vanishing result. Nevertheless they enter the spin sum rule on equal 
footing, which is satisfled in the model as we demonstrated. Furthermore, we have shown that also the subleading 
GPDs satisfy the polynomiality condition, and that the coefficients in front of the highest power in ^ in even moments 
of {H'" + H'^)(x, ^, t) and (E^ + E'^){x, t) arc equal to each other up to an opposite sign. Given the different technical 
handling in the model, the fulfillment of general requirements where both GPDs are involved is a strong check of the 
consistency of the approach. 

In the chiral quark-soliton model = 2 J'^ already at a low scale, i.e. the fractions of nucleon momentum and 
spin carried by quarks and antiquarks are the same. In QCD this relation becomes exact in the limit of a large 
renormalization scale. Actually this is not a specific prediction of the chiral quark-soliton model, but should hold 
in any model describing the nucleon consistently in terms of quark and antiquark degrees of freedom, i.e. without 
gluon degrees of freedom. For such models to be consistent M'^ = 2 J"^ = 1 must hold. The chiral quark-soliton 
model fulfills this requirement. In the framework of the xQSM the contribution of gluons to the nucleon momentum 
and angular momentum is suppressed with respect to the corresponding quark contributions by the small parameter 
characterizing the packing fraction of instantons in vacuum [40] . 

We then focused on the GPD (S" + £"^)(.x, ^, and computed numerically its forward limit with a constant 
(momentum-independent) constituent quark mass M using the Pauli-Villars regularization. We found that (i?" + 
E'^){x,0,Q) is negative at larger x and sizeable, and positive at small x, changing the sign at a; w 0.2, while (S" -|- 
E'^)(x, 0, 0) is always positive and concentrated towards smaller x. Noteworthy is that in the; model in the chiral limit 
{E^ + E'^){x, 0, 0) is strongly enhanced at small x with respect to (if" -|- H'^){x, 0, 0). It would be interesting to study 
in more detail the chiral mechanism responsible for that. 

On the basis of our results we estimated azimiithal asymmetries in hard exclusive production of vector mesons and 
pion pairs from a transversely polarized proton target. We observe that pion pair production is very sensitive to 
(£J" -|- E'^){x, ^, t). Data from experiments at HERMES or JLAB on this reaction would be well suited to distinguish 
among different models. 

Our work in particular the demonstration of the consistency of the approach - not only increases the faith into 
the applicability of the chiral quark-soliton model to the computation of GPDs, but also practically paves the way to 
further studies of GPDs in leading and subleading order of the large- A^c limit. 
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Note added in proof 

After this work was finished Ref. [60] appeared, where among others a fit of {E'> — E^){x,0,t) to the Pauli form 

factor was reported. The result of Ref. [60] for (_£" + E'^ — i?" — E'^){x, 0, 0) is qualitatively similar to the ansatz of 
Ref. [4] plotted in Fig. 4c and differs from our result which is negative at all x. 

APPENDIX A: MODEL SYMMETRIES 

The effective Hamiltonian HeS, Eq. (11), commutes with the parity operator II and the grand-spin operator 

K = S L -h t (Al) 

where S = 5757°7, L = f x p and T = denote the single-quark spin, angular momentum and isospin operators, re- 
spectively. Thus, the single-quark states are characterized by the eigenvalues of II, if^ and as \n) = ]i^„, tt, iv', K^). 
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Consider the unitary matrix Gs given in the standard representation of Dirac- and flavour SU(2) matrices by 
G5 = 7^7"''r^. It acts as G^^^G'^^ = (7'')^ and GsT^G^^ = — (r")"^ such that in coordinate space the Hamiltonian 
(11) is transformed as G^HG^^ = . This transformation acts as G5$„(x) = $*(x). The G5 transformation in 
the effective chiral theory corresponds to time reversal transformation accompanied by a flavour rotation [61]. Using 
the G5-transformation one flnds {m\T°'\j) = — (j|T°|m) and, e.g., 



(m|(l+7V)i^(p)G(X)i/(p)|j) = (j|(l-7V)i?(-p)G(X)F(-p)|m) 
(m|r^(l + 7V)i^(p)G(X)7J(p)|j) = -(j|r^(l - 7V)i?(-p)G(X)i^(-p)|m) . 



(A2) 



Another useful symmetry of the model is parity n|TO) = Pm\fn) with Pm = ±1 depending on the parity of the state. 
We obtain, e.g.. 



(m|r^(l+7V)i^(p)G(X)i?(p)|i) = p„p,(m|r^(l - 7V)^^(-p)G(-X)F(-p)|i) . 



(A3) 



APPENDIX B: POLYNOMIALITY 



For the discussion of moments the following expressions for HE{x,^,t) and EM{x,^,t) are more convenient 

1 



V m,oc< 



E 



-iz"Ei 



m,all 

JiOCC 



Em — Ej 



121 J 2tt 

X {j\ T^il + 7V) exp(-i0V/2) exp(iAX) exp(-i^V/2) |m) 

+tz'^ e-*"''^"(m|3(l +7°73) exp(-i2V/2) exp(iAX) exp(-izV/2) |m) 



v V m,oc< 



X (mlr'-lj) (j| (1 + 7"7') exp(-i2V/2) 



E 

m,all 

J,OCC 



exp(iAX) exp(-i2V/2) |m) 



(Bl) 



+-°E 



(m|T''(l + 7''7-') exp(-izV/2) 



^3ab 



exp(iAX) exp(-izV/2) |m) 



(B2) 



If we take out from the double sums the limiting case Em Ej and combine it with the respective single sums in 
(Bl, B2) then we arrive at the expressions in Eqs. (28, 29). 

1. HE{x,i,t) 

In order to demonstrate that ait = the moments of He{x, t) are in fact functions of ^ only, we have to continue 
analytically the expression in Eq. (Bl) to the unphysical point t = Q. For that note that 



exp(iAX) = (-'2gM^)'' |x|'.p,Jcosg^) 



(B3) 



ie=0 



where cos^ = X^/|X| and Pi^ are Legendre polynomials, cf. Ref. [17]. Taking moments of HE{x,^,t) in Eq. (Bl) 
and using (B3) we obtain 



M 



1 

(^,0) = jdx x'^-^HE{x,i,t) 



TT 



E 

.Li=0 



J=0 



m.all 

JiOCC 



Em — Ej 



Li 



7=0 



E 



(B4) 
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with 



rnj 

,E _ 



(mlr'-li) (jlr^Cl +7V) {fY |X|'eP;,(cos^) {f) 

m) 



m) 



^ (m|3(l + 7V) {fY |Xr^P;,(cose) {f) 



3\J 1 



f?\Li-J 



Using the Gs-symmetry, cf. Eq. (A2), we obtain 



^ (-l)^*0-|rV) {m\ t\1 - 7V) {fY |Xi''P,,(cos^) {fY^-' \j) , 
Bl ^ (-l)^^(m|3(l - 7%') {fY \±\''PiA'=os6) {fY''' \m) • 



(B5) 



(B6) 



Exploring these identities (for the first it is necessary to rename under the sum in (B4) m ^ j) we see that and 
B^ in Eq. (B4) are effectively given by (note that (7*^7^)^ = 1 for even N and 7O73 for odd N) 



<i = {m\A3) (il T^(7°7')''* {fY \±\''PiA^ose) {fY'-' \m) , 
B^ = (m|3(7V)''^ {fY \±\'^PiAoose) {fY'-' \m) . 



(B7) 



Prom parity transformations, cf. Eq. (A3), we conclude that A'^j ~ = if /e is odd. Thus we observe that only 
even powers in ^ contribute to the moments. In the next step we demonstrate that the infinite series in is actually 
a polynomial. 

For that let us first observe that the object r*" \j) (j Ir** (where summation over b is implied) is invariant under hedgehog 
rotations. Thus it is practically an irreducible spherical tensor operator of rank zero with respect to simultaneous 
isospin- and space-rotations [17]. So in both cases, A^j and Bj^, we deal with matrix elements of the type {m\0\ni), 
the only difference being the "replacement" T''|j)(j|r*' 3 which is irrelevant for our arguments. The matrix elements 
{m\0\m) are non-zero only if the operator O is rank zero, cf. [17] for details. O is a product of 7''7^ which have 
rank 1, and P;^ which is rank 1^, and rank zero operators. O is in general a reducible operator which, however, can 
be decomposed into a sum of irreducible tensor operators. For our purposes it is important to find the highest value 
for le in the Legendre polynomial which yields a non- vanishing result for (m|0|m). The operators p", which appear 
Xj-times in whatever ordering, and operator 7°7'^, for odd Li, can combine to an operator with maximally rank Li + 1 
for odd Li {Li for even Li). Thus le can have at most rank Li + 1 for odd Li {Li for even Li), otherwise the total 
operator O cannot reac;li the needed rank zero. Considering that Li < {L — 1) we see that the L*^ moment of 
is a polynomial with the highest power for even L {^^~^ for odd L). Thus we obtain 



Mi^ne,o) = 



MnNc 



E 

Xi=0 

L-2 



ii=0 



E ( J ) ^ E 

J=0 ^ ^ le=0 
le even 



Li 



J=0 



!e=2 
ie even 



(-^2a/.v) 



I, -2 



m.occ 
J, all 



m.all 

JiOCC 



Em. — Ei 



(B8) 



2. EM{x,i,i) 

In this case the relevant structure to be continued analytically to t = is, cf. Ref. [17], 

-z2CMw)'=-2 



|f exp(iAX) = 



Ze=2 



[P1,|X|'»P;^(C0S^)] . 



(B9) 



With this result we obtain for the moments of Em{x, ^,t) at zero t the result 

1 



M 



{L) 
M 



{^,Q) = jdxx'^-^EM{x,i,t) = ^^Y.- 



1^-2 



(^-i)E(\f)E(";)^E^^-^-^^" 

r n ^ ^ r— n ^ ^ m,occ 



Li jL^\ J 

1,4=0 ^ ^ 7=0 



Li=0 ^ ^ J=0 ^ ' ^ "•?,"=■= "".all 

J, all J,OCC 



(BIO) 
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with 



{m\T'\j) 01 (1 +7%') e'^^'lpl, Pi Acose)] (#)^^-^ |m) , 



= (m|r^l+7V)(r)-'e3«^[pl,|X|'»Pi,(cos^)](# 



f,3\Li-J 



m) . 



(Bll) 



Prom the Gs-symmetry, cf. Eq. (A2), we obtain 



'-mj 



(B12) 



Prom parity transformations, cf. Eq. (A3), it folfows that and B^ are zero if Ig is odd. Thus we observe again 
that only even powers in ^ contribute to the moments. 

In order to demonstrate that the infinite series in even powers of ^ is actually a polynomial, let us consider first the 
double sum contribution in Eq. (BIO). We deal with the matrix elements which are of the type (m|r"'|j)0|O''|m). 
The operator r" is an irreducible spherical tensor operator of rank 1 with respect to simultaneous isospin- and space- 
(hedgchog-) rotations [17]. Thus only certain transitions (™|T°|j) between the states |m) and \i) are non-zero. In 
order for the entire expression for A^^ to be non-zero the same transitions must appear also in the piece {j\0°'\m), 

thus the operator 0° must be rank 1, too. 

The operators r°, 7*^7^ and P;^, which in our context arc relevant for the rank counting, can combine to an 
operator with maximally rank Li + 2 for even Lj (Lj -|- 1 for odd Li). Thus 1^ can be at most Li -\- 2> for even Li 
{Li + 2 for odd Li), otherwise the total operator cannot reach the needed rank 1. However, 1^ must be even, thus 
le ^ {Li + 2) for even {Li + 1 for odd Li). Finally Li < {L — 1). Therefore the double sum contribution to the L*^ 
moment of M^^^ is a polynomial with the highest power Jq,. gyen L (^^^^ for odd L). 

B^ in Eq. (B12) is zero unless the operators (|X|, r°, the Li + 1 operators p", Pi^, and j'^'y^ if Li is odd) combine 
to a rank operator. One obtains the same restriction for the maximal value of le in term of Lj as above, however, 
in this case Li<{L- 2). Thus, we finally obtain (with A^,, B^ given in Eq. (B12)) 



E 

.Li=0 



L-l\^fLA 1 ^ (-/2ii\Lv) 



E ( J ) 2^-. E 

J=0 ^ ' ie = 2 



/,,-2 



-(^-i)E 



Li=0 



L-2 
Li 



le even 

Li +2 



m,all 

j , occ 



.1 



M 



Em. — Ei 



E ( ,7 ) oLi E 



{-t2^M, 



N 



\l.-2 



J / 2^i 

,7=0 ^ ^ 'e=2 



(B13) 



We know that (i?" +H'^){x, ^, t) satisfies polynomiality and that its even moments L are polynomials of degree L [17]. 
Here we observe that even moments L of Em{x.S,, t) are polynomials of degree (L — 2). Since Em{x, t) is the sum of 
{E" + £"^)(.T, and {H" + H'''){x, ^, f) this means that the moments of (£"" + _B'')(.x, ^, i) also satisfy polynomiality 
and that the coeSicients of the highest power of even moments satisfy the relation (5). 



APPENDIX C: PROOF OF THE SPIN SUM RULE 



The model expression for the second moment of EM{x,^,t) is given by (B13) for L = 2. By observing that 
\p1,X'^P2{cos9)] = iXl we obtain 



M^^\i,0) = -^ 



E 



{Ern + Ej){m\T'\j){j\{j^^^)e^-'Xl\ 

Em — Ej 



E 



2{m\T^\j){j\p^e^-bXl\m) 



Em — Ej 

Using the following identity, where we explore the equations of motion in the model, 

{j\ {Em+Ej)j°j'Xi\m) = {j\ {h, 7V ^i} |m) = {j\ (2p' X^ + ij^ j') |m) 

we obtain 



2-^ r^cc L^rn - Ej 

j,non 



{m\T'\j) {j\ 7%'75|m) + 2(m|r^|j) (j|e^'^''X V|m) 



(CI) 



(C2) 



(C3) 



19 

In Eq. (C3) wc identify the model expressions for the spin S*^ and angular momentum L*5 contribution of quarks and 
antiquarks to the spin of the nucleon 

j,non 

^'^ = E {j\^'''x^p'H , (C4) 

-' m.ooc 
j , non 

where denotes the isosinglet axial coupling constant. The operators in (C3) can be added according to Eq. (Al) 

^7°7^75 + e'^^'^X'f = S'^ + P = K^-f^ . (C5) 
Noting that in Eq. (C3) the matrix elements {j\K^\m) = Kj6jm vanish we obtain 

m1^)(C,0) ^ 2SQ + 2LQ = ^Y1 ^r4^(m|r3|j) {j\r^\m) = 1 . (C6) 

m.occ 



[1] D. Miillcr, D. Robaschik, B. Goycr, F. M. Dittos and J. Horcjsi, Fortsch. Phys. 42, 101 (1994). 
A. V. Radyushkin, Phys. Lett. B 385, 333 (1996); Phys. Rev. D 56, 5524 (1997). 
X. D. Ji, Phys. Rev. D 55, 7114 (1997). 

,]. C. Collins, L. Frankfurt and M. Strikman, Phys. Rev. D 56, 2982 (1997). 
[2] X. D. Ji, J. Phys. G 24, 1181 (1998). 
[3] A. V. Radyushkin, arXiv:hop-ph/0101225. 

[4] K. Goeke, M. V. Polyakov and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 47, 401 (2001). 

[5] M. Diehl, Phys. Rept. 388 (2003) 41. 

[6] X. D. Ji, Phys. Rev. Lett. 78, 610 (1997). 

[7] J. P. Ralston and B. Pire, Phys. Rev. D 66, 111501 (2002). 

M. Burkardt, Int. J. Mod. Phys. A 18, 173 (2003). 
[8] M. V. Polyakov, Phys. Lett. B 555, 57 (2003). 

[9] For an overview sec: M. Vanderhaeghen, Nucl. Phys. A 711, 109 (2002). 
[10] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 87, 182001 (2001). 
S. Stepanyan et al. [CLAS Collaboration], Phys. Rev. Lett. 87, 182002 (2001). 

C. Adloff et al. [HI Collaboration], Phys. Lett. B 517, 47 (2001). 

S. Chekanov et al. [ZEUS Collaboration], Phys. Lett. B 573 (2003) 46. 
[11] M. Vanderhaeghen, P. A. Guichon and M. Guidal, Phys. Rev. D 60, 094017 (1999). 
N. Kivel, M. V. Polyakov and M. Vanderhaeghen, Phys. Rev. D 63, 114014 (2001). 
A. V. Belitsky, D. Miiller, A. Kirchner and A. Schafer, Phys. Rev. D 64, 116002 (2001). 
V. A. Korotkov and W. D. Nowak, Eur. Phys. J. C 23, 455 (2002). 
A. V. Belitsky, D. MuUer and A. Kirchner, Nucl. Phys. B 629, 323 (2002). 

A. Freund and M. F. McDermott, Phys. Rev. D 65, 074008 (2002); Eur. Phys. J. C 23, 651 (2002). 
A. Freund, M. McDermott and M. Strikman, Phys. Rev. D 67, 036001 (2003). 
A. Kirchner and D. Mtiller, arXiv:hep-ph/0302007. 
[12] M. V. Polyakov and A. G. Shuvaev, arXiv:hep-ph/0207153. 

[13] P. V. Pobylitsa, Phys. Rev. D 65 (2002) 077504 and 114015; Phys. Rev. D 66 (2002) 094002; Phys. Rev. D 67 (2003) 

034009 and 094012; arXiv:hep-ph/0211160. 
[14] X. D. Ji, W. Melnitchouk and X. Song, Phys. Rev. D 56, 5511 (1997). 

[15] V. Y. Petrov, P. V. Pobylitsa, M. V. Polyakov, I. Bornig, K. Goeke and C. Weiss, Phys. Rev. D 57, 4325 (1998). 
[16] M. Penttinen, M. V. Polyakov and K. Goeke, Phys. Rev. D 62, 014024 (2000). 

[17] P. Schweitzer, S. Boffi and M. Radici, Phys. Rev. D 66, 114004 (2002); Nucl. Phys. A 711, 207 (2002). 

P. Schweitzer, M. Colh and S. Boffi, Phys. Rev. D 67 (2003) 114022 [arXiv:hep-ph/0303166]. 
[18] A. Mukherjee and M. Vanderhaeghen, Phys. Lett. B 542, 245 (2002); arXiv:hep-ph/0211386. 

[19] S. Boffi, B. Pasquini and M. Traini, Nucl. Phys. B 649 (2003) 243; arXiv:hep-ph/0311016. 

[20] A. Radyushkin, Armalen der Pliysik 13 (2004) 718 [arXiv:liep-ph/0410153]. 
[21] M. Vanderhaeghen, Annalen der Physik 13 (2004) 740. 

[22] D. I. Diakonov and V. Y. Petrov, JETP Lett. 43 (1986) 75 [Pisma Zh. Eksp. Teor. Fiz. 43 (1986) 57]. 
[23] D. I. Diakonov, V. Y. Petrov and P. V. Pobylitsa, Nucl. Phys. B 306, 809 (1988). 

D. I. Diakonov, V. Y. Petrov and M. Praszalowicz, Nucl. Phys. B 323 (1989) 53. 

[24] C. V. Christov, A. Z. Gorski, K. Goeke and P. V. Pobylitsa, Nucl. Phys. A 592 (1995) 513 [arXiv:hep-ph/9507256] . 
[25] C. V. Christov et al, Prog. Part. Nucl. Phys. 37, 91 (1996). 
[26] D. I. Diakonov et al., Nucl. Phys. B 480, 341 (1996). 



20 



[27] P. V. Pobylitsa and M. V. Polyakov, Phys. Lett. B 389, 350 (1996). 

[28] D. Diakonov, V. Y. Potrov, P. V. Pobylitsa, M. V. Polyakov and C. Weiss, Phys. Rev. D 56, 4069 (1997). 

[29] C. Weiss and K. Goekc, arXiv:hcp-ph/9712447. 

[30] P. V. Pobylitsa et al, Phys. Rev. D 59, 034024 (1999). 

[31] D. Diakonov, V. Y. Petrov, P. V. Pobylitsa, M. V. Polyakov and C. Weiss, Phys. Rev. D 58 (1998) 038502. 

[32] M. Wakamatsu and T. Kubota, Phys. Rev. D 60, 034020 (1999). 

[33] K. Goeke et al, Acta Phys. Polon. B 32, 1201 (2001). 

[34] P. Schweitzer et al, Phys. Rev. D 64, 034013 (2001). 

[35] M. V. Polyakov and C. Weiss, Phys. Rev. D 60, 114017 (1999). 

O. V. Teryaev, Phys. Lett. B 510, 125 (2001). 
[36] D. L Diakonov and M. L Eides, Sov. Phys. JETP Lett. 38, 433 (1983) [Pisma Zh. Eksp. Teor. Fiz. 38, 358 (1983)]. A. Dhar, 

R. Shankar and S. R. Wadia, Phys. Rev. D 31, 3256 (1985). 
[37] D. L Diakonov and V. Y. Petrov, Nucl. Phys. B 272, 457 (1986) 457. 
[38] D. L Diakonov and V. Y. Petrov, Nucl. Phys. B 245, 259 (1984). 
[39] E. Witten, Nucl. Phys. B 223, 433 (1983). 

[40] D. Diakonov, M. V. Polyakov and C. Weiss, Nucl. Phys. B 461 (1996) 539 [arXiv:hep-ph/9510232]. 
[41] J. Balla, M. V. Polyakov and C. Weiss, Nucl. Phys. B 510 (1998) 327 [arXiv:hep-ph/9707515]. 

B. Dressier and M. V. Polyakov, Phys. Rev. D 61 (2000) 097501 [arXiv:hep-ph/9912376] . 
[42] See Sec. 3 of Ref. [4]. Cf. also P. V. Pobylitsa and M. V. Polyakov, Phys. Rev. D 62 (2000) 097502. P. V. Pobylitsa, 

arXiv:hep-ph/0203268. A. V. Efrcmov, K. Goekc and P. V. Pobylitsa, Phys. Lett. B 488 (2000) 182. 
[43] X. D. Ji, J. Tang and P. Hoodbhoy, Phys. Rev. Lett. 76 (1996) 740 [arXiv:hep-ph/9510304]. 
[44] P. Hagler and A. Schafer, Phys. Lett. B 430 (1998) 179 [arXiv:hep-ph/9802362]. 

A. Harindranath and R. Kundu, Phys. Rev. D 59 (1999) 116013 [arXiv:hep-ph/9802406]. 
P. Hoodbhoy, X. D. Ji and W. Lu, Phys. Rev. D 59 (1999) 014013 [arXiv:hep-ph/9804337]. 

[45] M. Wakamatsu and T. Watabe, Phys. Rev. D62, 054009 (2000) [hep-ph/9912500]. 

[46] O. V. Teryaev, arXiv:hep-ph/9904376, and arXiv:hep-ph/9803403. 

[47] A. Freund, Eur. Phys. J. C 31 (2003) 203 [arXiv:hcp-ph/0212017]. 

[48] S. Kahana and G. Ripka, Nucl. Phys. A 429 (1984) 462. 

[49] M. Wakamatsu and H. Yoshiki, Nucl. Phys. A 524 (1991) 561. 

[50] A. Freund and L. M. Sehgal, Phys. Lett. B 341, 90 (1994). 

S. D. Bass, S. J. Brodsky and I. Schmidt, Phys. Lett. B 437, 417 (1998) [arXiv:hep-ph/9805316]. 
[51] T. Kubota, M. Wakamatsu and T. Watabe, Phys. Rev. D 60 (1999) 014016 [arXiv:hep-ph/9902329]. 
[52] M. Gliick, E. Rcya and A. Vogt, Eur. Phys. J. C 5 (1998) 461; Z. Phys. C 67 (1995) 433. 

[53] B. Lehmann-Dronkc, A. Schafer, M. V. Polyakov and K. Gooke, Phys. Rev. D 63, 114001 (2001) [arXiv:hep-ph/0012108]. 

B. Lehmann-Dronke, P. V. Pobylitsa, M. V. Polyakov, A. Schafer and K. Goeke, Phys. Lett. B 475 (2000) 147 [arXiv:hep- 
ph/9910310]. 

[54] A. V. Radyushkin, Phys. Rev. D 59, 014030 (1999) [arXiv:hep-ph/9805342]. 

[55] E. J. Brash, A. Kozlov, S. Li and G. M. Huber, Phys. Rev. C 65, 051001 (2002) [arXiv:hep-ex/01 11038]. 
[56] P. E. Posted, Phys. Rev. C 51, 409 (1995). 

[57] I. V. Musatov and A. V. Radyushkin, Phys. Rev. D 61, 074027 (2000) [arXiv:hep-ph/9905376]. 

[58] Wc used the QCD evolution program of QCDNUM available at http://www.nikhef.nl/ h24/qcdnum. 

[59] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 28 (2003) 455 [arXiv:hep-ph/0211080]. 

[60] M. Diehl, T. Feldmann, R. Jakob and P. KroU, arXiv:hep-ph/0408173. 

[61] P. V. Pobylitsa, arXiv:hep-ph/0212027. 



